Are Cluster Magnetic Fields Primordial ? 
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We present results of a detailed and fully non-linear numerical and analytical investigation of 
magnetic field evolution from the very earliest cosmic epochs to the present. We find that, under 
reasonable assumptions concerning the efficiency of a putative magnetogenesis era during cosmic 
phase transitions, surprisingly strong magnetic fields 10 -13 — 10 -11 G, on comparatively small scales 
100 pc - 10 kpc may survive to the present. Building on prior work on the evolution of magnetic 
fields during the course of gravitational collapse of a cluster, which indicates that pre-collapse fields 
of ~ 4 x 10 -12 G extant on small scales may suffice to produce clusters with acceptable Faraday 
rotation measures, we question the widely hold view that cluster magnetic fields may not be entirely 
of primordial origin. 



Magnetic fields exist throughout the observable Uni- 
verse. They exist in the interstellar medium, in galax- 
ies, and clusters of galaxies (for reviews cf. 0). The 
origin of galactic- and cluster- magnetic fields is still un- 
known. A plausible, though by far not convincingly es- 
tablished possibility is the generation of magnetic seed 
fields and their subsequent amplification via a galactic 
dynamo mechanism. Seed fields may be due to a vari- 
ety of processes (and with a variety of strengths), such 
as the Biermann battery within intergalactic shocks 0, 
stellar magnetic fields expelled in planetary nebulae, or 
during supernovae explosions, either into the intragalac- 
tic, or in the presence of galactic outflows into the inter- 
galactic medium , as well as due to quasar outflows of 
magnetized plasma Seed fields may also be of pri- 
mordial origin with a multitude of proposed scenarios. 
These include generation during first-order phase transi- 
tions (e.g. QCD or electroweak) , around cosmic defects, 
or during an inflationary epoch (with, nevertheless, ex- 
tremely small amplitudes) . For a review of proposed sce- 
narios we refer the reader to p|. 

The philosophy in prior studies of primordial magneto- 
genesis is often (but not always) as follows. After estab- 
lishing a battery mechanism (e.g. separation of charges 
and production of currents) and a "prescription" or es- 
timate for the final, non-linearly evolved magnetic field 
strength (e.g. equipartition of magnetic energy with tur- 
bulent flows), subsequent evolution is approximated by 
simply assuming frozcn-in magnetic field lines into the 
plasma. Though such an approximation may be appro- 
priate on the very largest scales, it should be clear, that 
this may not be the case on the "integral" or coherence 
scale of the field. Here, coupling of the magnetic fields to 
the gas induces non-linear cascades of energy in Fourier 
space. The characteristics of initially created magnetic 
field are thus vastly modified during cosmic evolution be- 
tween the epoch of magnetogenesis and the present. 

The final step in such studies is then often to deter- 
mine field strengths on some prescribed large scale (e.g. 
lOMpc) typically falling in the range 10~ 30 G < B < 
10~ 20 G, inferring that this may act as seed for a suffi- 
ciently efficient dynamo to produce galactic- and cluster- 



magnetic fields of order 1CP 6 G. This is observed in neg- 
ligence of the fact that much stronger fields on smaller 
scales, result not only from a variety of astrophysical 
seeds, but from these very same primordial scenarios. 
Considering the likelihood of a magnetized plasma in the 
early Universe, it seems important to be able to make 
predictions on such final remnant fields surviving a mag- 
netogenesis scenario to the present epoch, particularly 
also since such fields may fill voids of galaxies and may 
potentially be observable by upcoming technology [||. 

The purpose of this letter is twofold. We have at- 
tempted to develop a coherent picture of gross features 
of non-linear, cosmic MHD evolution of primordial fields, 
including all relevant dissipative processes, such as vis- 
cosity due to diffusing- or free-streaming- neutrinos and 
photons, as well as ambipolar damping. A subset of the 
results of our numerical and analytical analysis is pre- 
sented here, whereas details are presented elsewhere 0. 
Our study allows us, for the first time, to make predic- 
tions for magnetic field energy and coherence length at 
the present epoch, for broad ranges of initial magnetic 
configurations, parametrized by spectral index, initial 
helicity, initial energy, as well as era of magnetogene- 
sis. Second, drawing on a prior numerical MHD study of 
the gravitational collapse of a cluster of galaxies , we 
challenge the often-cited conclusion that cluster magnetic 
fields may not be entirely of primordial origin. Rather, 
we stress that it seems not clear at the moment if such 
fields on comparatively small scales may not, after all, 
produce cluster magnetic fields as widely observed. 

In passing we note that there exists a number of ana- 
lytical HEfllUHIIlQEI , and numerical H[l7|, 
studies on the evolution of non-helical and helical pri- 
mordial magnetic fields, which, nevertheless, for one or 
the other reason either remain inutile in predicting final 
field properties, or do so only for a specific scenario. The 
inutility of results is related to facts such as, the adop- 
tion of an evolutionary model not supported by numerical 
simulations, an inadequate treatment of viscosity due to 
photons, or simply, the analysis being linear in nature 
or being performed in Minkowski space and not properly 
transferred to the expanding Universe. 
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FIG. 1: Decay of magnetic energy in different damping 
regimes and for different inital conditions as observed in our 
numerical simulations. The t~ 12 line shows the theoretical 
damping law in the turbulent regime for n — 3, h g — 0. 

The generation of primordial magnetic fields in magne- 
togenesis scenarios is generally believed to occur during 
well-defined periods (e.g. QCD-transition). Subsequent 
evolution of these magnetic fields should therefore be well 
approximated by a "free decay" without any further in- 
put of kinetic or magnetic energy, i.e. as freely decay- 
ing MHD. The exceedingly large Prandtl numbers in the 
early Universe allow one to neglect dissipative effects due 
to finite conductivity. Nevertheless, as already stressed in 
Ref. ^3 1 dissipative effects arising from the "imperfect- 
ness" of the fluid due to neutrino- and photon- diffusion 
and free-streaming play an important role in early MHD 
evolution. One may further show that for field strength 
as considered in this letter the assumption of fluid incom- 
pressibility is appropriate. 

To verify theoretical expectations we have performed 
numerical simulations of incompressible, freely decay- 
ing, ideal, but viscous MHD. These simulations are per- 
formed with the help of a modified version of the code 
ZEUS-3D [Hll9| in a non-expanding (Minkowski) back- 
ground and on 128 3 to 512 3 grids. Modifications lie in 
the inclusion of fluid viscosities [20|. It is shown else- 
where 0, Il6j| that conformal- or near-conformal- invari- 
ance of the MHD equations allow for the interpretation of 
results obtained with Minkowski metric to results appli- 
cable for a Friedman-Robertson- Walker metric. Results 
of such simulations, in particular the decay of magnetic 
energy i? ma g with time, for a variety of physical regimes 
are shown in Fig. ^ 

We have found that results of our simulations may be 
understood in a comparatively simple manner. In partic- 
ular, non-linear MHD processing of the initial spectrum 
at epoch with Hubble constant H(T) occurs for all scales 
I, which obey 

v(l)/l>H(T), (1) 

irrespective of the Reynolds number R e of the flow. Here 
v(l) may be written as the Alfven velocity vaQ) = 
-^(0 / V '^(g + p) when turbulence holds, R e > 1, and 



as v(l) = va(1) R e (v = va,1) for viscous MHD (R e < 1). 
This holds equally during the photon diffusion (Z 7 <C I) 
regime R e — vl/ij and the photon free-streaming (l^ 3> I) 
regime R e = v/al, where rj, ev, and l 7 are photon shear 
viscosity, drag coefficient, and mean free path, respec- 
tively. Defining L(T), the integral- or coherence- scale, 
as the scale where equality applies in Eq. Q), one finds 
that L(T) is in fact, the scale containing most of the en- 
ergy of the flow. This is due to a non-linear and rapid 
cascade developing on all scales I < L, with energy in 
fluid eddies transported down to the dissipation scale ^diss 
and transferred to heat. Since the resultant small-scale 
spectrum is red and we assume the initial as yet unpro- 
cessed large-scale (I > L) spectrum to be blue, L(T), as 
the smallest unprocessed scale, remains as the magnetic 
coherence scale of the field at epoch with temperature 
T. Magnetic energy is then approximately given as the 
initial energy on scales L(T). 

For fields which are maximally helical, i.e. TL = 
H m ax ~ {B 2 (l) I) « B 2 (L) L, Eq. (TJ may still be used to 
obtain the coherence scale of the field. Nevertheless, due 
to helicity density TL = (1/V) J d 3 ir A • B (where A is the 
vector potential and V is the integration volume) being 
an ideal invariant in the early Universe, a direct cascade 
of energy from large scales to small scales is accompa- 
nied by an inverse cascade of energy from small scales to 
large scales. That is, whereas in the absence of helicity 
the large-scale field remains unprocessed, in the maxi- 
mally helical case large-scale fields undergo growth even 
on scales I > L{T). Surprisingly, we find that during this 
process of large length scale magnetic field amplification 
the initial spectral index n is conserved [13. The decay of 
magnetic energy, is thus described by the requirement of 
conservation of helicity, in conjunction with an increase 
of L{T) described via Eq. Due to a vast increase of 
L(T) in the early Universe, even initially sub-maximally 
helical fields, i.e. TL g = h g TL max . g , with h g < 1, ulti- 
mately reach a maximal helical configuration. Here, and 
throughout, an index 11 g" denotes properties at the mag- 
netogenesis epoch. Parameterizing the initial (comoving) 
magnetic energy spectrum by qb ~ QBg{l/ L g )~ n , where 
L g is the initial (magnetogenesis) coherence scale obeying 
Eq. and gs g is approximate initial magnetic energy, 
one finds that fields have reached a maximally helical 
state when L(T) has grown beyond 

This picture may be employed to derive damping ex- 
ponents, i.e. .©mag ~ £~ 7 , and compare to those in- 
ferred from numerical simulations (cf. Fig.^). Whereas 
the comparison is quite favorable in the viscous regime, 
turbulent decay is observed somewhat slower than pre- 
dicted. For example, for non-helical, turbulent MHD 
with a n — 3 spectrum we predict 7 = 1.2 whereas 
the best fit to the numerical simulation yields 7 = 1.05. 
Nevertheless, we argue that this trend, seen also by oth- 
ers ^3 El- must be due to limited numerical resolu- 
tion. Its explanation, i.e. a putative additional increase 
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FIG. 2: Evolution of comoving coherence length for different 
initial magnetic field configurations with r fl = 0.083 and n — 3 



(solid lines from top to bottom: h g = 1, h g 
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The mean free paths of neutrinos and photons are labeled by 
l v and l~f, and Ih is the Hubble length. 



of relaxation time beyond l/v with continuing evolution, 
requires the existence of additional dimensionful quan- 
tities (of dimension velocity or length) associated with 
the flow. As we may find none, other than i^iss (which is 
mostly widely separated from L in the early Universe, but 
not in the simulation) or Lhoxi the size of the simulation 
box, we believe the effect to be numerical in nature. We 
have noted, that spectra at late times show a peak region 
AL quite spread, and are likely only marginally resolved 
by the simulations. In any case, larger numerical simu- 
lations are required to address this effect, with resulting 
predictions for the surviving magnetic fields, given in this 
letter, being on the conservative side. 

We have undertaken the in practice straightforward 
but arduous effort to assemble these results and, under 
the inclusion of all appropriate dissipation terms, and for 
quite general initial conditions, followed the growth of 
magnetic coherence length and energy density from the 
very earliest times to the present jjj. Fig. [3 shows exam- 
ples for the growth of L(T) for a number of scenarios of 
magnetogenesis at the QCD phase transition. The evo- 
lution is observed as an alternation between turbulent 
MHD and viscous MHD. "Viscosity" here is early on due 
to neutrinos, some time before recombination due to pho- 
tons, and after recombination due to hydrogen-ion scat- 
tering and hydrogen-hydrogen scattering. Particularly 
notable are phases where the growth of L(T) is halted 
completely. This occurs either at epochs before recom- 
bination in the viscous regime with diffusing photons or 
neutrinos, as well in part of the regime when those parti- 
cles are free-streaming or at epochs after recombination, 
due to the peculiar redshifting of Eq. and/or the ef- 
fects of hydrogen diffusion and ambipolar drag. Note, 
however, that the growth of L(T) and concomitant de- 
crease of B(T) during the late phases of viscous MHD 
with free-streaming photons (neutrinos) may be faster 



than the growth of L(T) during turbulent MHD. Note 
also that initial conditions leading to relatively strong 
magnetic fields at recombination result in a rapid in- 
crease of L(T) at T rcc ps 0.3 eV, whereas for weaker fields 
B < 1CP 13 G a similar jump occurs at reionization. Note 
that effects due to structure formation are not taken into 
account here (see below, however). 

We give here the final (pre-structure formation) results 
on the coherence scale L(T ) and field amplitude B(T ), 
where To is present CMB temperature, derived by em- 
ploying Eq. QJ, as we U as retaining the initial field energy 
due to all scales I > L in the submaximally helical case, 
and conserving helicity density in the maximally helical 
case. Complete results on L and B for all eras may be 
found in [Jj. Fields which remain still submaximally he- 
lical (i.e. L(T ) < L™ ax ) [22| at the present epoch have 
for coherence length and field strength 



B{T ) k 1.65 x 1(T 6 G x n/{n+2) 



X (0.O83) 
L(T ) « L gc x- 2 /^ ( 



V2 / T„ 



100 MeV 
100 MeV 



-n/(n+2) 



2/(n+2) 
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where x = 2.30 x 10 9 is a small factor and r g = 

{qb / Sr^)g is a convenient measure of magnetic energy 
density in terms of radiation entropy density s r = 
(4/3) g g (7r 2 /30) at the magnetogenesis epoch assumed 
to occur at temperature T g . Note that, somewhat opti- 
mistically r g = 0.083, when magnetogenesis results in 
magnetic energy density equivalent to the photon energy 
density shortly after a QCD-transition with g g »s 10.75. 
The comoving coherence length L gc at the magnetogen- 
esis epoch is given by 



l„* 0.45 pc v^yy 



1/2 



T„ 



100 MeV 



(5) 



This yields, for example, for r g = 0.083, T g = 100 MeV, 
n = 3 to the appreciable field strength of B c « 1.1 x 
10~ n G. If we were to apply the simulation observed 
instead of theoretically predicted decay exponent 7 the 
result would increase approximately to B ~ 5 x 10~ n G. 
On the other hand, fields which have reached a maximally 
helical (i.e. L(Tq) > L™ ax ) state at present have 



B c w 4.69 x 1Q- 12 G y 
L c m 550 pc y ^/n 



(6) 
(7) 



with 



r 9 V/ 2 ( h g 
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where h g is the fractional helicity of the maximal one 



at the generation epoch. 
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FIG. 3: Final magnetic field strengths in the r g — h g param- 
eter space for T g = 100 MeV (solid lines) and T g = 100 GeV 
(dotted lines), where we assumed a spectral index of n = 3. 

Though on small scales, it is seen that surprisingly 
strong fields may survive an early Universe magnetoge- 
nesis epoch to the present (cf. Fig. I3J. This is inter- 
esting in light of recent simulations on the formation 
of clusters of galaxies from slightly overdense and pre- 
magnetized regions via the gravitational instability 
The authors find, that fields of strength B c « 4 x 10~ 12 G 
(corresponding to their quoted B « 10~ 9 G at simulation 



starting redshift z = 15) yield to clusters whose Faraday 
rotation measures are essentially indistinguishable from 
those observed in real clusters |23|. Furthermore, the 
authors arrive at the intriguing conclusion, that this re- 
sult is virtually independent of whether a homogeneous 
field is assumed initially, or a field whose energy contri- 
bution is dominated by fluctuations on the very smallest 
scales in their simulations |2^|. It is not clear if such an 
erasure of memory of initial conditions, possibly related 
to an interplay between the development of shear flows 
and small-scale turbulence during the course of gravita- 
tional collapse, pertains if initial field coherence lengths 
in the cluster simulations are reduced by a further factor 
~ 100 (due to the comparison between typical coher- 
ence lengths in Eqs. J3J, Q and the spatial resolution, 
~ 100 kpc comoving, of the simulations). Nevertheless, 
if so, cosmological magnetic fields generated during early 
eras, either of moderate magnetic helicity, or generated 
fairly late, could account for present-day observed cluster 
magnetic fields, and as such in the absence of any further 
dynamo amplification. This holds true even if as recently 
claimed such fields have very blue sprectra |2^| . We con- 
clude that this interesting possibility seems to deserve 
further investigation. 
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